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Abstract 

We present some recent progresses on Heun functions, gathering results from classi¬ 
cal analysis up to elliptic functions. We describe Picard’s generalization of Floquet’s 
theory for differential equations with doubly periodic coefficients and give the de¬ 
tailed forms of the level one Heun functions in terms of Jacobi theta functions. The 
finite-gap solutions give an interesting alternative integral representation which, at 
level one, is shown to be equivalent to their elliptic form. 


1 Introduction 

Heun functions [9] are defined as a natural generalization of the hypergeometric function, 
to be the solutions of the Fuchsian differential equation 

r d^F Jj d ee \ dF (g + afdlCuj) ^ ^ 

/ dw‘^ \w 1 — w 1 — kJw) dw w{l — w){l — kJw) ’ 

with four regular singularities 0, 1, 1/A;^, oo. In [3] the fourth singularity is a = 1/fc^ and 
the auxiliary parameter s is taken as q = —sjUJ. We will follow the notation 

s; a,/9, 7 , (5; tc), fc^G[0,1]. (2) 

In [3] one can hnd a lot of information on these functions and in the remaining chapters 
of [21] are gathered many results on their confluent limits. Due to limitation of space 
the former section misses several important items, particularly the relation between Heun 
functions and elliptic function theory. These last years this held has received new inter¬ 
esting developments from workers in condensed matter physics [5] and integrable systems 
[15]. Interesting accounts have been given by Smirnov [23] and by Takemura [24]. 






Our aim is to gather, for the community dedicated to the study of special functions 
and orthogonal polynomials, these new progresses relating Heun and elliptic functions. 

As a preliminary step we will give one further motivation for the study of Heun func¬ 
tions, coming from orthogonal polynomials and birth and death processes. 

2 Prom orthogonal polynomials to Heun functions 

Let us consider the three terms recurrence [1] 

xPn bji—i Pji—X ~\~ CliYiPn bji u ^ 1. (3) 


We will denote by Pn and Qn two linearly independent solutions of this recurrence with 
initial conditions 


Pq{x) = 1, Pi{x) = 


X — ao 


^0 


Qo{x) = 0, Qi{x) = 


1 


(4) 


The corresponding Jacobi matrix is 

/ ao bo 0 0 0 ■ ■ ■ \ 

bo ai bi 0 0 ■ ■ • 

0 bi 02 &2 0 


V 


(5) 


J 


If the bn > 0 and G M the Pn (resp. the Qn) will be orthogonal with respect to a 
positive probabilistic measure (resp 


Pm{x) Pn{x) dllji^x) bnin'j 


QmQQ QnQQ dlf)^ Errin') 


( 6 ) 


with the moments 


Sn= X"^ dlp{x), U > 1, <§0 = 1 


(7) 


When considering applications to birth and death processes with killing [14], it is 
sufficient to consider the special case where 

T ^^n T Jn; -^n/^n+l i ^ ^ 0, 

AqAi . . . A„_i 

TTo = 1, TTn = -, U > 1. 

/^l/^2 ' ' ' 

It is convenient to introduce new polynomials Fn{x) by 

(-ly 


Pn{x) = 


Fn{x), 


n >0, 


( 8 ) 





The positivity constraints 


^ O5 hn+1 ^0 ^ 0 ? 

will be assumed. From (3) we deduce 

^rF'yj Hn+iF^j^x (A„ T /i„ T T 

F_i(a;) = 0 , Fo(a;) = 1 . 

The orthogonality relation ( 6 ) becomes 


(9) 


( 10 ) 


j Fm{x) Fn{x) d'lp{x) = 7r„(5, 


( 11 ) 


The particular case where Xn-, p-n are quadratic is connected with Heun functions in 
the following way. Let us take 


Xn = k'^in + a){n + P), = n(n + 7 — 1 ), 7 „ = /c'^(5n 

a > 0, /5 > 0, 7 > 0, /c^ G (0, 1 ), k'"^ = 1 — k'^, 

and consider the generating function 

F{x] w) = Fn{k‘^, X] a, f3, 7 , 6 ) 10 "^ = Hn{k‘^, s; a, /3,^, S;w), s = x — a^k"^. 

n>0 


( 12 ) 


(13) 


Routine computations show that F(a;; w) is the solution, analytic in a neighbourhood of 
to = 0 of Heun’s differential equation ( 1 ). Quite remarkably, the most general quadratic 
birth and death process with linear killing produces, for its generating function, the most 
general Fuchsian equation with four regular singularities. 

Let us conclude with the proof of: 


Proposition 1 The Hamburger (hence the Stieltjes) moment problem corresponding to 
the polynomials Fn{x) with recurrence coefficients ( 12 ) is determinate. 

Proof : 

A theorem by Hamburger [l][p. 84] states that if the series Cn = Pn(0)^ is divergent then 
the Hamburger moment problem is determinate, i. e. the measure ip dehned in ( 6 ) is 
unique. From relation ( 8 ) we can write 

c«+l _ /in+l Fn+l{0f _ (n + l)(n + 7 ) Fn+l(0)^ 

Cn Xn Pn(0)2 {u + a) {u +/ 3 ) Fn{ 0 )‘^ 

Now the generating function F{x-,w), introduced in (13), is analytic for |t(;| < 1: hence 
its radius of convergence is one. It follows that, for n ^ 00 the previous ratio has for limit 
and the series {cn} diverges. □ 

To have a clear view of the problems encountered in the construction of solutions of 
Heun’s equation, we will introduce some terminology: we will call “generic” solutions the 
solutions valid for arbitrary values of the auxiliary parameter s and “non-generic” solutions 


^In [21] this function is called “local” Heun. Notice that we have for normalization Hn{- • •; 0) = 1. 



those which require particular values of s. Let us begin with some results on non-generic 
solutions. 


Non-generic solutions 


1 Special hypergeometric cases 


The reason for the restriction G (0,1) is that in the two limiting cases = 0 (notice 
that the positivity conditions (9) require k"^ > 0) and k"^ = 1 the four singular points 
reduce to three and therefore Heun functions degenerate into hypergeometric functions. 
Indeed we have: 

1. For k^ = 0 the parameters {a, j3) become irrelevant and we have: 


{ iLn( 0 ,s; (a,/3), 7 , 5; w) = 
r± = p ± V/FTs, 

This case, which would correspond to a —cx), is 


2^1 


r+, r_ 
7 




P = 2 ( 9 ' + '^- !)• 
missing in [3]. 


(14) 


2. For = 1 we have the relation: 

Hn{l, s] a, 13, 7 , 3; w) = (1 - wY 2 F 1 


r + a,r + {3 
7 


; w 


(15) 


= p - Vp2 - 0,3 - s, P = 2(7 - a - /3)' 


This case is considered in [3] but only with the extra constraint s = —aj3, i. e. for r = 0. 


2 The “trivial” solution 

This solution corresponds to the special values s = a(3 = 0. In this case Heun’s differential 
operator is factorized as 

Dyj -1- M F = 0, 

with obvious L and M. This factorization leads to 

= ^ U = = (16) 

Starting from these solutions we can use the change of function 

F^F: F{w) =w^{l-wY{l-k‘^wyF{w) (17) 

which transforms, for special values of the parameters p, a, r, a Heun function into another 
Heun function, up to changes in the parameters [3][p.l8]. In such a way one generates 7 
more Heun functions starting from (16). For all these solutions Heun’s differential operator 
remains obviously factorized and, as shown in [22], this happens only for these cases. 



3 Derivatives of Heun functions 


In general the derivative of a Heun function cannot be expressed in terms of some Heun 
function with different parameters. However, as shown in [11], this happens in 4 cases: 

Hn'{k‘^, 0; a, /?, 7 , 5; w) = — Hn{k‘^, s'; a + 2, /? + 2 ,7 + 2, 5 + 1; w), 

7 + 2 (18) 

-s' = 7 + 5 + 1 +(7 + 6 + l)k‘^. 

Hn'lk"^, —a/3k‘^; a, f3, 7 , 5; w) = — —k‘^(l — w) Hn(k‘^, s'; a + 2, /? + 2 ,7 + 1, 5 + 2; w), 

7 + 1 

-s' = 7 + 5 + 1 +(7 + 6 + aj3)k‘^. 


Hn'ik'^, —tt/d; a, f3, 7 , 5; w) = (1 — k'^w) Hn{k‘^, s'; a + 2, /? + 2 ,7 + 1, 5 + 1; w), 

7 + 1 

-s' = 7 + 5 + Q;/ 5+(7 + e + l)k‘^. 

( 20 ) 

Hn '(/c^, s; 0, /?, 7 , 5; w) = - Hn{k‘^, s'; 2, /? + 1 ,7 + 1, 5 + 1; w), 

7 + 1 (21) 

s' = s — 7 — 5 — (7 + e)k‘^. 

Notice that, from (13), the functions Hn{- ■ ■ ;tc) are analytic around tc = 0. 


(19) 


4 Reduction to hypergeometric functions 

For non-generic solutions it was realized some time ago the possibility for Heun functions 
to reduce to hypergeometric functions. Some relations, using Weierstrass elliptic functions 
are given in [13]. Later on Kuiken [16] has observed some reductions to hypergeometric 
functions of some particular rational variable R{w). This may happen only for polynomial 
R{w), with the following list : 

Hn (—1, 0; 2a, 25, 2c — 1,1 + a + 5 — c; tc) = 2+1 ^ 

Hn (1/2, —2a6; 2a, 25, c, 1 + 2(a + 6 — c); tc) = 2+1 ^ ;'?^(2 —w)), ( 22 ) 


Hn (2, —4a5; 2a, 25, c, c; w) = 2+1 


a, 5 


; Aw{l — w) 


Notice that in all these cases R{w) is a second degree polynomial with 3 free parameters. 
These results have been completed recently by Maier [17], who discovered new cases ^ with 
cubic and quartic polynomials R{w) and 2 free parameters: 


Hn (1/4, —9a5/4; 3a, 35,1/2, 2(a + 5); tc) = 2+1 ^ ~ > 

Hn (1/2, —8a5; 4a, 45, a + 5 + 1/2, 2(a + 5); tc) = 

^•^'(0 + 6+1/2 i M 2 -!»)(!- + 


( 23 ) 


^We have omitted the cases involving complex values of 5^. 



These results are not the whole story, since Heun functions may reduce to the product of 
some function f{w) by some hypergeometric function with variable R{w), not necessary 
polynomial. A result of this kind was obtained in [ 12 ] by Joyce : 

Hn (1/4, -1/8; 1/2,1/2,1,1/2; w) = 2^1 

with 

R{w) = ^ (^2 — tcVd — w — (2 — tc)Vl — , tc e [0,1]. 

Let us turn ourselves to the generic solutions of Heun’s equation. 


Generic solutions 


1 The 192 solutions of Heun’s equation 

As Heun himself observed [9] there is a set of 24 substitutions of the variable w which 
produce a transformation of Heun’s equation into another Heun’s equation with differ¬ 
ent parameters. This leads to a complete list of 192 solutions. This list has been fully 
worked out recently by Maier in [18]. We will just quote the generalizations of the Euler 
transformation of the hypergeometric function: 

s; a, P; 7 , S, e; w) = 

(25) 

= (1 — —s — 7(5 — 1 ); a — J -|- 1, /5 — J -|- 1 ; 7 , 2 — J; tc), 

and of the Pfaff transformation 
Hn{k‘^, s] a, /?; 7 , 6, e; w) = 

= (1 — Hn{—kP/k' —k‘^/k'‘^{s + 0 : 7 ); q;,q; — ( 5 -|- 1 , 7 ,q; — /5-|-l;t(;/(t(; — 1)). 

(26) 


2 An integral transform 


An integral transform was given in [26]: 


HniJP, s] a, P, 7 , J; w) = 

r(7) 


r(a)r(7 - a) Jo 


t°‘ ^( 1 —t)''' " ^ Hn{k'‘^, s;'y, P, a, 6 +'f — a;wt) dt, 


(27) 


valid for Re 7 > Rea > 0 and w E C\ [1, cx)[. 

Let us use orthogonal polynomials to prove this relation [27]. To this end we dehne 
the monic polynomials Mn by 


Mo = Fo = 1, Mn{P;x) = =n\{j)nFn{P;x), n>l, 



where P denotes the set of parameters a, (3, 7 , S. These monic polynomials satisfy the 
recurrence relation; 

(■^n T hn T ''fn ^n+1 T -^n—ihn ^n—li ^ ^ 0, 

= 0 , Mo{x) = 1 . 

If we dehne the new set of parameters 

P' = (/c^, a' = 'y, (3' = (3, 'y' = a, 6' = 6 + 'f — a), x' = s + a' 
it is easy to check the invariance relation Mn{P'',x') = Mn{P',x). This induces 

fo;) 

x = s + aPk‘^, Fn(P]x)= Fn(P';x'), x'= s + a'(3'k'^. (28) 

(7)n 

For Re 7 > Re q; > 0 we can write 

_ _ ^( 7 ) _ f .n+a—1/-| j.\ 7 —a—1 

( 7 ). “r(a)r(7-a) io ^ ^ 

then multiplying both sides of this relation by (28) and summing for n > 0 gives (27) for 
|t(;| < 1. Analytic continuation extends it to tc G C\ [1, cxd[. 

Remarks : 


1. This is not an integral equation, since the parameters of the Heun function are 
changed in the transformation. Notice that several integral equations are known [3]. 

2. For 0 we recover Bateman’s integral relation 


F(A)F(c-A) io 


^A-1(1 _ a,b ^29) 


valid for Rec > ReA > 0. Notice that now A is a free parameter, and this enables one to 
deduce Euler’s integral representation. This does not work for Heun functions because the 
parameter a is not free. 

3. For k"^ = 1 we get the relation 


{I-WY2F, 


r + a,r + j3 


■w = 


r(7) r 

F(a)F(7-a) Jq 


p-\i - {1 - wty 2 F 1 


r + 7 ,r + 


; wt J dt^ 


^ 

r = p- _aj3-s, p = -{^-a-f3), f = p-- s, p =-a +j3). 


This relation does not appear in the extensive list of hypergeometric integrals given in [7], 
so it could be new. 



3 Carlitz solutions 


In his analysis of some orthogonal polynomials of Stieltjes, Carlitz [4] discovered the follow¬ 
ing remarkable result: the linearly independent Heun functions with parameters (/c^, s 7 ^ 
0 ; 0 , 1 / 2 , 1 / 2 , 1 / 2 ) are given by 


exp ( ± 2i ^/sz{w)), 


du 

ziw) = / — , 

Jo 2 a / m (1 — m )(1 — k?‘u) 


(31) 


To check most simply this result we note that by the inversion theorem of elliptic functions 
we can write 


z{w^ k'^) = 


du 


lo 2 a/m(1 — m)(1 — k’^u) 

This conformal transformation maps the singular points as 


w = sn^(; 2 , k"^). 


w 


0 

i 


1 


00 


K 


K + iK' iK' 


where K = K{k‘^) and K' = K{k'‘^) are the complete elliptic integrals of hrst kind. Using 
2 : as new variable and setting F{w) = y{z), Heun’s equation becomes 


A 

dz"^ 


(27 -1) 


cnzdnz 
sn 2 


(2S - 1) 


sn 2 dn 2 
cn 2 


(2e- 1)A;2 


sn 2 cn 2 
dn 2 


dy 

dz 


-|-4(s 4 - a(3 k"^ sn ‘^z)y = 0. 


(32) 


At the very symmetric point where a(3 = 0 and 7 = <5 = e = 1 / 2 , it reduces to 




which proves (31). 

Using then the transformation (17) we can generate a full set of 8 generic solutions (see 
[26]) where the detailed list is given). So we realize that Heun’s equation, in the elliptic 
functions setting, lies in the held of differential equations with doubly periodic coefficients. 
It happens that Floquet’s theory for differential equations with periodic coefficients does 
generalize to the case of doubly periodic coefficients and was derived by Picard. This will 
be discussed in the next subsection. 


4 Elliptic functions of second kind and Picard’s the¬ 
orem 

Since this theorem is not easily available in the standard textbooks [2],[28], we will present 
some background material. Let us just recall the dehnition and a few basic results from 
elliptic function theory : 



Definition 1 A function $( 2 :) is elliptic if it is meromorphic and it has two linearly 
independent periods 


^{z+ 2uj) = ^{z), ^{z+ 2uj') = ^{z). 

These two periods define therefore a non-degenerate period parallelogram (usually ui = 
K, u' = iK'). 

We will need also 

Proposition 2 An elliptic function: 

1. Has as many poles as zeroes in a period parallelogram. 

2. Can he written 

, Hiz-bA-'-Hiz-bn) 

— A— ---—- ai + ■ ■ ■ + On — hi + ■ ■ ■ + bn, (33) 

H[z — ai) ■ ■ ■ H[z — On) 

where H{z) is one of Jacobi’s theta function defined as 

H(z) = 2 sin ((2n + 1)|^) , , = e-w. 

n>0 


3. Must have at least a pole of multiplicity 2 in a period parallelogram otherwise it 
is a constant. 

4. With poles of multiplicity n at z = a can be expanded as: 

H'iz) 

^{z) = '^aog{z - a) + aig\z - a) ^ - ang^'^\z-a), 

where the sum is extended to all the poles. 


To state Picard’s theorem we need the dehnition 


Definition 2 A function <h(; 2 ) is a an elliptic function of second kind (or a function with 
constant multipliers) if it is meromorphic and that there exists two non-vanishing constants 
/i and p' such that 

^z-J 2uj) = ^z) ^{z + 2uj') = e^^' ^z), 

where 2uj and 2uj’ are two linearly independent periods. The constants pi and pi' are called 
the multipliers. 

Obviously for = 1 we recover elliptic functions. We will have to discuss 

separately the two cases: 

1 . Generic multipliers for which up,' — u'p 7 ^ 0, 

2 . Special multipliers for which up' — u'p = 0. 



4.1 Theorems for generic multipliers 


Let us consider the function 


f{z) = e 


XzH{z-a) 


( 35 ) 


H{z) ’ 

where a is not homologous to zero. One can check that it is elliptic of second kind with 
generic multipliers 

/i = 2XK, fi' = i (^XK' — • 

Let us prove: 


Proposition 3 An elliptic function of second kind with generic multipliers: 

1 . Admits as many zeroes as poles in a period parallelogram. 

2. Must have at least a simple pole in a period parallelogram, otherwise it vanishes. 

3. Can he expanded, using f{z) defined in (35), as 

= X] -a) + cif'{z - a) H-h Cnf‘^'^\z - a), (36) 

where the sum extends to all poles z = a of multiplicity n of <h. 

Proof : 

Notice that for any elliptic function of second kind *h(;2) with given generic multipliers 
(/i, fj,') it is always possible to hnd A and a in (35) such that / has the same multipliers. 
It follows that ^{z)/f{z) is elliptic. Using (33) we get the most general structure 

^ / X . X H{z — hi) ■ ■ ■ Hiz — hn) 

<h(^) = A f{z) — -r-—- ai + ■ ■ ■ + On = hi + - \-bn, 

H[z — ai) ■ ■ ■ H[z — On) 

which proves the first assertion. 

For the second assertion, let it be supposed that <h has no pole: it cannot have any 
zero and therefore *h(;2) reduces to Ae^^, but in this case the multipliers are not generic, 
contradicting our hypothesis, hence A = 0. 

For the third assertion, let us consider a pole z = a of multiplicity n. Near toz = awe 
have the Laurent expansion 

f^’^Uz) = - -+ holomorphic, 

[z — 

so we can hnd coefficients {c^, k = 0 .. .n} such that 

$( 2 :) - '^{cof{z -a) + cif'{z - a) H-h Cnf'"\z - a)) 

has no poles in a period parallelogram, so it must vanish. □ 



4.2 Theorems for special multipliers 

This time let us define 

M ( 37 ) 

One can check that it is elliptic of second kind with special multipliers. We have now 

Proposition 4 Any elliptic function of second kind with special multipliers can he ex¬ 
panded as 

^{z) = C e^^ + ^ {cof{z - a) + Ciffz - a) H- (z - a)) , (38) 

with the constraint 

e ^“(co + Cl A + ■ ■ ■ A"^) = 0, (39) 

where the sums extend to all poles of ^{z). 

Proof: 

In this case u p,' — u' fj, = 0. So we can find a valne of A such that F{z) = ^{z)e~^^ is 
elliptic. Let us consider a pole 2 : = a of order n of <h. We have for Laurent’s series 

<h( 2 :) = — --—H-h (-iTj — . + holomorphic. 

z — a [z — a)^ [z — a)"-+^ 

hence we can write 

F{z) = -+ ■ ■ ■ + (-I)"" . + holomorphic, 

z — a [z — ay [z — 

where the new residue is 

Oq = e '^^(co + Aci + ■ ■ ■ + \^Cn). 

Using the expansion theorem for elliptic fnnctions (34) we have 

<h(^)e = C + e oq — a) + oi ~ + ■ ■ ■ 

- OnD^ ~ • 

Inserting e^^ into the right-hand side, and expanding the derivatives according to Leibnitz 
rnle gives (38). The constraint (39) comes from the fact that the sum of the residues in a 
period parallelogram vanishes for an elliptic fnnction. □ 


4.3 Picard’s theorem 

Now we can state Picard’s theorem [19] 



Proposition 5 Let us consider the differential equation 


d^F , ^d^-^F 

with doubly periodic coefficients 

ai{z + 2u!) = ai{z), ai(z + 2u') = ai(z), l = l,2,...,n, 

Any meromorphic solution F{z) is a linear combination of elliptic functions of the second 
kind. 


Proof : 

To shorten, and in view of application to Heun’s case, we will consider a differential 
equation of second order: 

d‘^F dF 

^+p{z)^ + q{z)F = 0, (40) 

dz^ dz 

with doubly periodic coefficients p and q. Let us consider a solution F{z) which is not 
elliptic of second kind: then the ratio F{z+2uj)/F{z) cannot be a constant. The periodicity 
of the coefficients implies that the functions F{z + 2uj) and F{z + 4a;) are also solutions, 
so we must have a linear relation of the form 


F{z + 4a;) = AF{z) + BF{z + 2a;). 


Let us now consider the non-vanishing function (j){z) = F{z + 2u) + pF{z). If we take for 
p a root of -|- Bp — = 0 it is easy to check that (f){z + 2uj) = {B + p)(f){z). 

A similar argument works with respect ot the period 2a;'. So we have proved that we 
can construct a first solution (j){z) which is elliptic of the second kind. 

We have now to prove that a second linearly independent and meromorphic solution, 
which can be written as 

ij{z) = (j){z) J G{z)dz, P' = P, 

is also elliptic of second kind. 

The Wronskian 

fj (j)' — fj' (j) = C e“^, C* 7^ 0, 

and the meromorphy of 0 and tf imply that e~^ is meromorphic. Since p is doubly periodic, 
e~^ is elliptic of second kind. So G is a meromorphic elliptic function of second kind. 

Let it be supposed first that G has generic multipliers. Using the expansion theorem 
(36) we can write 


= X] Co/(^-a)H- \-Cn f^'^\z - a) , 


f{z) = e 


Xz H{z-a) 

H{z) ■ 


where the sum includes all the poles of multiplicity n of G{z). Since the integral has to be 
meromorphic all the coefficients Cq must vanish and we get 


Q{z)= jG{z)dz = '^ cif{z-a)^ - ^ Cn ^\z - a) , 



and from (36) it follows that Q{z) and 'ip{z) are elliptic of second kind. 

Let us now consider the case where G has special multipliers. The expansion theorem 
(38) gives 

G(z) = Ce^^ + J2[cofiz-a) + --- + g f^^\z - a)] , /(z) = 

with the constraint 

^(co + cih + --- + c„h”)e-'^“ = 0. (41) 

If h 7 ^ 0 the meromorphy of G requires that all coefficients Cq vanish, so we can write 

g{z) = ^ ^ [ci /(^ - a) H-h {z - a) . 

This relation and the constraint (41), with coefficients cq all vanishing, implies that g{z) 
is elliptic of second kind with special multipliers, hence ip is again elliptic of second kind. 
If h = 0 we have 

g{z) = Gz + '^^cif{z - a) -\ - \-Cn - a) , ^ cq = 0. 

This implies 

g{z + 2uj) = g(z) + D, G{z + 2uj') = g(z) + d' 

and so, if (Ai, A 2 ) are the multipliers of 0 we can write 

ip{z + 2uj) = Ai iIj{z) + XiD (j){z), il){z + 2u') = A 2 ip{z) + X 2 D' (j){z), 

so we can subtract from ip a suitable term linear in cp for which ip will be elliptic of second 
kind. □ 


5 The meromorphic solutions 

So let us look for the necessary conditions on the parameters to get meromorphic solutions 
in the variable ; 2 . The computation of the exponents at the singularities of (32) is quite 
simple and gives 

2: = 0 (sn2:)^^i (0,7 — 1) 

z = K {cnzy^^ (0,5 — 1) 

(42) 

z = K + iK' {dnzfP^ (0,e-l) 

z = iK' {snz)~‘^P^ (a,/3) 

So the necessary conditions for meromorphy are 


7 = - - mi. 


5=1- m 2 . 


e = 1 - m 3 . 


a = -^{mo + M), 0 = A(mo - M + 1), 


M = mi + m 2 + m 3 , 
iV = mo + M, 


( 43 ) 



with the vector N = (mo, mi, m 2 , m 3 ) G That these conditions are also sufficient is 
proved in [ 8 ]. From Picard’s Theorem the solutions will be elliptic functions of the second 
kind. The differential equation becomes 


(Py cnzdnz sn;2dn2; ,nsn2;cn2:\ dy 

^ + 2 -mi-+ m2-+ m3 k^— - / 

dz^ \ sn^: cn;^ an z J dz 

+ (4s + N{N — 2mo — 1) k'^ sn ‘^z)y = 0, 


(44) 


and for N = (n, 0, 0, 0) we are back to Lame’s equation. 

It is possible, extracting from y suitable factors of w, to relate the negative and positive 
values of the parameters rrii as summarized in [23] [p. 296], so from now on we will consider 
TV e NT 

It is interesting to get rid of the derivative in (44) by the change of function 


y{z) = isnz)^^ [cnz)^^ {dnzT^ Y(z) 


d^Y 

dz"^ 


(V(z)-A)Y 


(45) 


with 


+ "^ 2(^2 + 1)^ + m3(m3 + + mo(mo + 1) fc^sn^;^ 

A = As + (mi + m 2 )^ + k‘^{mi + 


(46) 


This equation has been considered by Darboux [ 6 ]. In this short article (3 pages) he 
claims: 


• that the product of two solutions is a polynomial, which we denote by 

of degree N in w = sn'^z and of degree g in a = As. The knowledge of this polynomial is of 
paramount importance as we will see later. Let us quote Darboux: “Une fois le polynome 
T determine, I’integration s’acheve, comme on sait, sans aucune difficulte.” 

• that for half-integer values of mi, m 2 , m 3 this equation can be integrated (for 
arbitrary A i. e. for what we call generic solutions). But neither detailed proofs nor the 
explicit forms of the solutions (even for the simpler case of integer rrii) were given. 

We would like to point out that even for Lame’s equation with half-integer values of n 
only non-generic solutions, due to Halphen and Brioschi, are known, so that the claim of 
Darboux seems questionable. 

To come back to Lame equation, let us mention that its solutions are meromorphic for 
integer n, and their general form, due to Halphen and Hermite, is given for instance in 
[28] [p. 570-575]. However, since one has to solve a set of n linear equations, only for low 
values of n everything can be made explicit. 

Interestingly enough, this differential equation has also appeared in the quite different 
field of integrable models, particularly KdV equation. Then relation (45) can be inter¬ 
preted as a Schrodinger eigenvalue problem, with eigenvalue A (equivalent to the auxiliary 
parameter s) and potential V{z) given by (46). 

The solutions corresponding to mi & Z are called “hnite-gap” solutions and sometimes 
V(z) is called a Treibich-Verdier potential [25], after their work on the subject. The name 
“hnite-gap” refers to the appearance of a hnite number of energy bands in the Bloch 
spectrum, a phenomenon discovered a long time ago by luce [ 10 ]. 



6 Elliptic level one 

solutions 


Proposition 6 

The level one (N 

= 1) solutions, with a 

= 4s, are given by: 

parameter 

solution y{z) 

constraint 

'F((T, w) 

mo = 1 

.zic.)H{z-uj) 

Q{z) 

dn = a — /c2 

a + k'^w — 1 — 

mi = 1 

zZ{uj) 0(^ ~ 

0W 

dn = a -I- 1 

aw + 1 

m2 = 1 

.zM0i(^-^) 

dn = a -f- 1 — 

A;2 a{l — w) + 1 — k"^ 

0W 

^zZ(co)Hl{z-u) 

ms = 1 

dn = a 

(j(l — k‘^w) — 1 + k' 

e{z) 


(47) 


9' 

with Z{z) = 

Proof: 

We will give the detailed proof for mi = 1, the other cases being analyzed similarly. So 
we start with 

(fy cn z dn z dy 


dz"^ 


sn 2 ) dz 


+ ay = d, 


(48) 


and we look for a solution y = Taking derivatives we get 




y" fy'^" 


y' H' , 0', , y « 

- = + 7=(t: 


+ /c^ sn^2:. 


yj sn‘^{z-p) 


Let us dehne the auxiliary function 


X{z) = 


- 2 


cnzdnz y' 1 

sn^: y SYi‘^{z — p) 


+ sn "^z + a, 


so that proving that y vanishes identically is equivalent to proving that y is indeed a 
solution of (48). 

For 2 ) = —iK' one can check that the poles in y cancel automatically. Imposing that 


the pole at 


0 is absent gives p 


— (p). The absence of the pole at 2 : 
H 


p is easily 


checked. So we know that y is bounded in a period parallelogram and since it is doubly 
periodic, it is bounded in all the complex plane: by Liouville theorem it is a constant. 
Imposing that y vanishes for z = K hxes the value of p by the equation sn^p = —1/a. 
Then we switch to the new parameter uj = p + iK' and we express the solution in terms 
of uj using the transformation theory of theta functions. 

We can now determine the spectral polynomial d/{a-,z) dehned, up to an overall con¬ 
stant factor, by the product of the two solutions of Heun’s equation. Let us do it first for 
mo = 1. The product of the two solutions is computed, using the transformation theory of 
theta functions, to the identity 


H{z — uj)H{z + oj) 
02(;2)02(o^)02(O) 


; 2 2 ; 2 2 
K sn 2 : — K sn o;. 


(49) 



which is a hrst degree polynomial with respect to tc = sn^;^ and to a = 1 + — 

k'^sn'^u. □ 


Remarks: 

1. Notice that a is uniquely dehned up to congruence. 

2. All the solutions have genus g = 1- 

3. The other linearly solution is obtained, for generic s 7 ^ 0 by the change 


-2:. 


4. The solution for mo = 1 is due to Hermite [28] [p.573] and the solution for m 3 = 1 
is due to Picard [20]. 

5. The hrst two formulas for 4/((T;tc) agree with the results given in [23]. 

6 . For special values of a, or equivalently of a, we may fail to get two linearly indepen¬ 
dent solutions. In this case factoring out the solution y given previously gives the second 
solution via a quadrature. Let us give some examples for the solution with mo = 1: 


a = 1 + k‘^ 


a = 1 


a = k"^ 


Asnz + Bsnz{ ^{z) + ^-^z 
H K 


4 r, i'Hi, . E-k'^K 

Acnz + B cnz | — [z) H-—- 2 : 

Hi K 


0 ' ^ 

Adnz + Bdnz ( -z^iz) + —z 


(50) 


01 


K 


7. The general structure of the elliptic solutions, for arbitrary level, are given in [ 8 ]. 
They suffer from the same defect as the general solution of Lame’s equation: they lead to 
really explicit expressions only for low values of the level N. 


7 Finite-gap solutions 

In [23], the following facts were used: 

1. Let us call 4/((T = 4s; tc) the product of two solutions of Heun’s equation and let 
us dehne p (resp. q) as the coefficient of the derivative (resp. of the function) in (1). Then 
T must be a solution of the third order differential equation 

T'" -h 3pT" -h {p' + 2 / -h dg)'!'' -h 2(g' -h 2pg)'F = 0. (51) 

For the meromorphic solutions of Heun’s equations, with parameters 

{ 7=i-mi, 6 = I - 1712, e=|-m3, M = mi-hm2 + m3, 

Q; = -i(mo + M), /3 = l{mo - M + 1), N = mo + M, 

the product is a polynomial of degree N with respect to the variable w: 

'Fg_jv(a, w) = ao{a) + ai{a) -h av(o'), 


(52) 



and a polynomial of degree g in the parameter a: 


= ho{w)a^+ hi{w)a^ - Vhg{w). (53) 

The leading term is given by a hrst order differential eqnation dne to the fact that a 
appears linearly only in the coefficients of T' and T of eqnation (51). We have taken for 
normalization 

b^^w) = w^^{i - wr^{i - k^wr^, 
as can be checked from (47). 

2. Then, following a method dne to Lindemann and Stieltjes [28] [p. 420] adapted to 
Heun’s case, one looks for a solntion of the form 

Inserting this ansatz into Henn’s equation gives on the one hand 

N' 

TV 

(recall that all the rrii are positive) and on the other hand: 


N{w) = 


w™i(l -w)™ 2 (i - k‘^w)^^ 
^/w{l — w)(l — k'^w) 



2 TT" - T' 2 + + 4gT2 


(54) 


showing that is of degree 25 f + 1 in a. The fact that z/ is a constant is easily checked by 
differentiating relation (54) and using the differential equation (51). 

So we conclude that Heun’s functions (the so-called “hnite-gap” solutions), with the pa¬ 
rameters given above, have for integral representation: 


'hg,Ar(A;tc) exp 



l)m 2 (l 



(55) 


For the level one solutions already considered the polynomial z/^((t) is given by 

mo = 1 {a — l)((j — k‘^){o' — 1 — /c^) 

mi = 1 a{a + 1)((T -|- k"^) 

(56) 

m2 = 1 a((T — /c^)((j + 1 — k'^) 

m3 = 1 a{a — l)((j -1-1- k"^) 

A partial list of the hnite-gap solutions, up to level N = 5, has been given in [23]. 

These results show that there do exist integral representations for Heun functions, but 
they are not “cheap”. Furthermore, as we will show now on a particular example, these 
integral representations are just a different dressing of the elliptic solutions: they are a 
kind of “algebraization” of the elliptic solutions, but are exactly the same analytic objects. 



8 Finite-gap versus elliptic solutions 


We will show, for the level one hnite-gap integral representation with mi = 1, that it 
does coincide with its corresponding elliptic solution given by (47). The other cases can 
be relateded using completely similar arguments. For this identihcation it is sufficient to 
consider w G [0,1]. We start from the data 

a = —sn= (t((T + 1 )((T +/c^), ^((t; ta) = ate + 1. 


So the factor appearing in the exponential of relation (55) is 

1 I" k'^ snu enu dnu w dw 

2 J {1 + a w)^/w{l — w){l — k'^ w) 

Under the change of variable w = sn‘^z, with G [0, K], it becomes an elliptic integral of 
the third kind, computed in [28] [p. 523]: 


± 


k^sn cu cn a; dn a; sn dz 

1 — sn ‘^ui sn ‘^z 


= ±- ( In 


Q{z — u) 
&{z + on) 


+ z Z{u) 


So, keeping the plus sign, we get for the exponential term in (55) 


.ZH ejz-u;) 

^eiz + uY 

And we need to multiply this by the square root of 

T oc 1 — sn sn = sn {k^ sn ^(o; — iK') — k'^ sn ‘^z'j , 

and upon use of the identity (49) we are hrst left with 

H{-uj + z + iK’) H{-uj -z + iK') 


4/ oc 


02(z) 


and after use of the transformation theory for the theta functions we end up with 

Q{z + ci;) Q{z — u) 


4/ oc 


02(z) 


Gathering (57) and (58) we conclude that the hnite-gap solution (55) reduces to 

0(z) ’ 

which is nothing that the elliptic solution given in (47). 


(57) 


(58) 


9 Conclusion 

The hnite-gap solutions or their elliptic counterparts solve Heun’s equation for all cases 
where these solutions are meromorphic functions of the variable z. It is interesting to notice 
that these progresses have been evolving in relation with integrability considerations. It 
is quite difficult to say what new ideas will require the generalization of these results to 
cover the non-meromorphic solutions. These much more difficult problems are left for the 
future. 
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